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Abstract 

The chiral phase transition induced by a charged scalar field is investigated numerically 
in a lattice fermion-gauge-scalar model with U(l) gauge symmetry, proposed recently as a 
model for dynamical fermion mass generation. For very strong gauge coupling the transition 
is of second order and its scaling properties are very similar to those of the Nambu-Jona- 
Lasinio model. However, in the vicinity of the tricritical point at somewhat weaker coupling, 
where the transition changes the order, the scaling behavior is different. Therefore it is 
worthwhile to investigate the continuum limit of the model at this point. 
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1 Introduction 



Scalar fields are mostly introduced into the quantum field theoretical models in order to 
trigger some symmetry breaking already on the tree level, as in the conventional Higgs 
mechanism. However, fundamental scalar fields can also play an important role in a dynam- 
ical symmetry breaking occurring only beyond the perturbative expansion. An example are 
strongly coupled Yukawa models, which can exhibit spontaneous chiral symmetry breaking 
(S^SB) even if the bare scalar potential in the action does not have the form of the classical 
Mexican hat 0-g. 

In this paper we continue the investigation of another model of this type, the lattice 
chiral symmetric xU4>4 model, considered in this context in ref. ||. The model consists of 
staggered fermion field Xi strongly coupled gauge field U and scalar field in four dimensions 
(4D). The gauge symmetry is compact U(l), U G U(l), and the model is confining at strong 
coupling. Both matter fields have unit charge and are, consequently, confined. In one phase, 
the Nambu phase, S^SB is generated dynamically by the vectorlike interaction between the 
gauge field U and the fermion field Xi making the gauge invariant condensate (xx) nonzero, 
in analogy to the S^SB in QCD. The Goldstone bosons are composed of x and x- 

The role of the fundamental charged scalar field is twofold. First, it shields the Un- 
charge of the fundamental fermion x, so that a composite neutral physical fermion state of 
the form F = <$x can exist asymptotically in spite of the confinement of the U(l)-charge. 
In the Nambu phase its mass mp is generated dynamically. This mechanism of fermion mass 
generation has been called shielded gauge mechanism in ref. || . 

Second, the scalar field has the tendency to suppress the S%SB and for very strong 
gauge coupling induces a new second order phase transition, at which the chiral symmetry 
is restored and amp, the fermion mass in lattice units, continuously approaches zero. In 
the scaling region of this phase transition in the Nambu phase, the model would describe a 
massive fermion in the continuum if it was nonperturbatively renormalizable. 

We initiate the investigation of this crucial question by means of numerical simulations of 
the xU<p4: model with dynamical fermions. The chiral phase transition is studied by varying 
suitably the gauge coupling (3 = 1/g 2 , the hopping parameter k and the bare fermion mass 
arriQ. Previous studies of the present |5|-[13| and similar ||14H19[1 models have indicated that 
the chiral phase transition, at which amp vanishes, is smooth in some interval of strong 
gauge coupling, including {3 = 0. At some value /3e of (3 this line changes the order at the 
tricritical point E and continues as a first order transition line for larger (3. We confirm 
these earlier results in high statistics simulations and localize with good accuracy the phase 
transition line and less precisely the point E — 0.64) on it. We investigate the scaling 
behavior of several observables in the vicinity of the transition. 

At {3 = 0, the xU<p4 model is identical to the 4D Nambu- Jona-Lasinio (NJL) model 
on the lattice with the same global chiral symmetry. That model has been investigated 
recently in much detail [20| and its nonrenormalizability has been confirmed. We compare 
data obtained for (3 > on lattices of various sizes and at different values of the bare fermion 
mass airtQ with those from the NJL model, using analogous analysis methods. We find that 
in the interval < (3 < (3e — e, e< 0.1, the behavior of the xU<p4 model is very similar to 
that of the NJL model. 

In the close vicinity of the tricritical point E the behavior changes significantly, however. 
The analysis methods, that work well at smaller /3-values, fail to describe the finite size 
effects and the amo-dependence of the data. Furthermore, a new correlation length, which 
is of the order of the lattice spacing or less at smaller /3-values, gets large here, indicating 
the presence of a new state in the spectrum of the model. This scalar bound state S - 



has no counterpart in the NJL model. 

From these observations we tentatively conclude that there are two possibilities how to 
approach the continuum limit of the xU<p4 model from within the Nambu phase in such 
a way that the fermion mass amp scales, and a massive fermion F is thus found in the 
continuum limit: 

1. At some point on the chiral phase transition line at (3 < (3e (perhaps even at negative 
(3). Here some generalization of the NJL model is obtained, possibly similar to those 
discussed in ref . GI[ . Either it is a nonrenormalizable theory like the NJL model itself, 



or a Yukawa-like model. 

2. At the tricritical point E. The spectrum would contain a massive fermion F, a massive 
scalar S and the obligatory Goldstone bosons tt. It could well be that the model is 
renormalizable and that the shielded gauge mechanism works here, as conjectured in 
ref. ||. Then the question would be, whether the resulting theory is Yukawa-like or 
whether some new interesting universality class exists here. 

Until now we have been able to localize the tricritical point E and to analyze the properties 
of the chiral phase transition outside the scaling region of this point. We found that the 
same methods of analysis do not work in the vicinity of E. The understanding of the point E 
itself requires substantially more data and new analytic tools for their analysis and is beyond 
the scope of this work. 

The outline of the paper is as follows: In the next section we define the lattice xU^a model 
and describe its phase structure at strong gauge coupling. Here we collect the quantitative 
results on the position of the phase transition lines (see fig. |T|). In sects. 3 and 4 we then 
describe in some detail the observed properties of the first order part (/3 E < (3 < /?t — 0.9) 
and second order part (/3 < /3e), respectively, of the chiral phase transition line. Our current 
preliminary experience with the vicinity of the point E is presented in sect. 5. Sect. 6 
contains summary and outlook. In the appendix we describe another tool used in this paper 
in addition to the hybrid Monte Carlo method: the microcanonical fermionic average method 



developed in refs. P2,E3fl and extended to the xU(pi model [24 . 



2 Lattice xU(f>A model at strong coupling 

2.1 The xU(J)a model 

The xU4>4 model is defined on the euclidean hypercubic lattice in 4D as follows: The stag- 
gered fermion field x of charge one leads to the global U(l) chiral symmetry of the model in 
the chiral limit, i.e. when the bare fermion mass mo vanishes. The gauge field link variables 
U are elements of the compact gauge group U(l). The complex scalar field <fi of charge one 
satisfies the constraint |0|=1. 
The action is 

S X U4> — S x + Su + , (2.1) 



where 



S x = l^y^VuxXx U Xitt Xx+ft - Ul_^Xx-n + arriQ^XxXx , (2.2) 



X /i=l 



S v = P t 1 - M^p}] , (2-3) 



4 

Here Up is the plaquette product of link variables U x ^ and rj^ = (— l) XlH The 
hopping parameter k vanishes (is infinite) when the squared bare scalar mass is positive 
(negative) infinite. The bare fermion mass mo is introduced for technical reasons, and the 
model is meant in the limit m = 0. 



2.2 Equivalence of the \U<\) model to the four fermion theory in 
the strong coupling limit 

At /3 = 0, the xU<p model in d dimensions can be rewritten exactly as a lattice four fermion 
model ||. In the path integral of the model ( [2.1| ) the scalar and gauge fields can be integrated 
out. This results in a purely fermionic model with the action 



#4f = - ^ Y i G XxXxXx+nXx+f, ~ IV»X {XxXx+fi ~ Xx+^Xx)] + ^ Y X * Xx ■ ( 2 - 5 ) 
x fJ.=l x 

Here 

J v Ju — JdUe^WU = h(2K), 

r = r(K) = — and (z.o) 
Jl Jt = JdUe 2KRe ^ = j Q ( 2/c ). 

The fermion field has been rescaled by y/r. The parameter r is an analytic function of k 
increasing monotonically from r(0) =0 to r(oo) = 1. 

The action (|2.5|) obviously describes a lattice version of the NJL model. The four fermion 
coupling parameter G is related to k via r: 

1 -r 2 

G = —-- . (2.7) 

From (|2.6|) one sees that G is decreasing monotonically with increasing k; G = oo at k = 0, 
and G = at k = oo. 



2.3 Phase diagram of the xU4>4 model 

The phase diagram of the x^04 model with U(l) gauge symmetry obtained in this work is 
shown in fig. [I], and the positions of the phase transitions are listed in table [3]. The used 
methods are described in the next sections. The phase diagram is consistent with the earlier 



results of unquenched simulations fllQ j|12l , |13l| . Qualitatively it agrees also with the quenched 
simulations R 0-0. The major difference is that the points E, T and C lie at smaller f3 
values than in the quenched case. 

The chiral phase transition of the NJL model at /3 — (point N), seen in the chiral 
condensate and the fermion mass amp extends to nonzero values of (3. As pointed out in 
ref. [0], this can be derived by means of a convergent expansion around (3 = 0. Therefore, 
some properties of the NJL model have been expected to persist also at small nonzero (3. 
We find this is so until (3 ~ 0.55. 

The model exhibits the Higgs phase transition (line ET and its nearly horizontal contin- 
uation to larger (3) seen in bosonic observables like plaquette and link energies. The striking 
feature is that, within the numerical precision, the chiral phase transition joins the Higgs 
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Figure 1: The phase diagram of the xU4>& model defined in eq. ( |2.1| ) at mo = 0. Three 
phases, the Coulomb, Higgs and Nambu phases, are found for k > and (3 > or (5 < 0. 
The emphasized points are: 

N: critical point of the NJL model, which is a special case of the xU4>a model at /? = 0, 
E: critical Endpoint of the Higgs phase transition line, 
T: Triple point, 

C: phase transition from the Confinement (at strong gauge coupling) to 

the Coulomb phase (at weak gauge coupling) in the model without the scalar field. 

The dashed line corresponds to a 2 nd order phase transition, full lines to 1 st order transi- 
tions. The dynamical fermion mass generation takes place in the Nambu phase. A detailed 
discussion of this phase diagram is given in ref. ||. 

The circles represent extrapolations of pseudo critical points to infinite volume and chiral 
limit. Also shown are the results obtained by means of the micro canonical fermionic average 
method on a 6 4 4attice (crosses) in the chiral limit. 
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0.64 
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Table 1: Estimated position of the phase transitions in the chiral limit and on infinite 
volume. 



phase transition line ET at the tricritical point E, forming a smooth line NET. The NE and 
ET parts are of second and first order, respectively. To our knowledge the interweaving of 
the chiral and Higgs phase transition is not understood theoretically. 

We concentrate on the Nambu phase, which is the area below the NET line. Here the 
chiral condensate (xx) is nonvanishing. Both x an d 4> fields are confined, in analogy to the 
quark confinement in QCD. The mass amp of the fermion state F = <$x is nonzero, and 
thus the dynamical mass generation occurs. 

The other phases as well as various limit cases of the phase diagram are discussed in 
much detail in ref. ||. It is helpful to keep in mind that in the quenched approximation 
the phase diagram corresponds to that of the U(l) Higgs model (the present model without 
fermion field), which has been thoroughly investigated |25]-2"8|. Even in this approximation, 



the chiral phase transition is observable in the chiral condensate, which is a nonlocal function 
of bosonic fields. 



2.4 Definitions of observables 

For a localization of the phase transition lines and a determination of the particle spectrum 
the following observables are used: 

The normalized plaquette and link energies, defined as 



Ep = ^?E Re ^p}> (2-8) 

El = w S M4£4,„<W , (2.9) 



where V = L 3 T is the lattice volume. These observables are useful for the localization of the 



Higgs phase transition. It is convenient |28| to use the perpendicular component of these 
energies, E± defined as 

E± — E L cos 9 + E P sin 9, (2.10) 

9(/3) being the slope of the chiral phase transition line at point (3. 

For the localization of the chiral phase transition line we measure the chiral condensate 
(xx) with the stochastic estimator method. 



Further we introduce the quantity 

9 ivi i \ 

(2.11) 



mo=0 

where Aj are the positive eigenvalues of the massless fermion matrix. This quantity equals 
the chiral susceptibility in the chiral symmetric phase. Of course, in the phase with S^SB 
it differs from the susceptibility, actually it diverges as 0(V). However, Q can be rather 
easily calculated on small lattices by means of the microcanonical fermionic average method, 
explained in the appendix, and used for an estimate of the position of the pseudocritical point 
at mo = 0. 

To calculate the mass am F of the physical fermion we consider the gauge invariant 
fermionic field 

F x = <f>tx*, Fx = <Px% B , (2-12) 

and determine numerically the corresponding fermion propagator in momentum space 

G FAB ( P ) = l^Y< ei ^ A)XF *Fv e ~ i{v+ * B)v ) ■ ( 2 - 13 ) 



The tta-, as well as the V ^ below, are defined in refs. p9[pO| . We look for the projections 

Gl( Pt ) = i±ti \G F (6, Pt )TA , G^(p t ) = ±txG F (0,p t ) , (2.14) 

l. J u=t 



which we fit by the formulas 



Gti.Pt) = — — f ^ (2-15) 

sm p t + {am F ) z 



Z F am F 
sm 2 p t + (am F ) 2 

Both fits and also those in configuration space give consistent results. The presented results 
are determined by means of eq. ( 2.15 ). 

Further we consider the fermion-antifermion composite states, the "mesons" . The time- 
slice operators for the mesons are given by 

O tt (*)=E s S^ar,«, (2-17) 

X 

with the sign factors s^f and other details given in ref. 0. We measure their correlation 
functions using point sources. 

We are also interested in the scalar and the vector bosons, which are present also in the 
gauge-scalar model without fermions. The corresponding operators ]26|] are 

x I i=l J 



^(0 = ^Yl lm {AtU(m,i<f> S+ lt} , « = 1,2,3. (2.19) 

X 

The masses of the scalar boson arris and the vector boson amy are computed from the 
corresponding correlation functions in momentum space. We note that the states S and V 
are expected to mix with the mesons carrying the same quantum numbers. In the QCD 
language these would be the a and u mesons, respectively. 
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Figure 2: Thermal cycles at (3 = 0.70 and (3 = 0.80. Both E± and (xx) have hystereses at 
the same place. 

3 ET line 

The ET line to the right of the point E (fig. [I]) is a first order phase transition line, and 
a continuum limit cannot be constructed here. We have investigated it mainly in order to 
verify the coincidence of the chiral and Higgs phase transitions. A better understanding of 
the ET line is also helpful for the localization and investigation of the tricritical point E. 

As known from the studies of the U(l) Higgs model the ET line (Higgs transition) 

is in the quenched approximation of first order in bosonic variables E^, arris and amy- In 
this approximation also the fermionic observables (xx) and amp show discontinuities at the 
same values of the coupling parameters as the bosonic variables do. We have found similar 
discontinuities in the full model (i.e. with dynamical fermions) too. This is demonstrated for 
(3 = 0.70 and [3 = 0.80 on the 8 4 lattice for am = 0.04 in fig. 0. Here E± and (xx) develop 
hystereses in thermal cycles in the k direction with low statistics per point (100 trajectories). 

The positions of the hystereses in bosonic and fermionic observables, E± and (xx) re ~ 
spectively, apparently coincide. This is so also for other values of am and on a 6 4 lattice. 
The hystereses shift slightly to lower K-values as am Q is decreased (by Ak < 0.02 between 
arriQ = 0.06 and 0.02). Their position is quite independent of the lattice size, as usual for 
sufficiently strong first order phase transitions. These two observations allow us to estimate 
the positions of the chiral phase transition in the am = limit (by linear extrapolation) 
and infinite volume (by neglecting a volume dependence). The results are indicated by two 
circles on the ET line in fig. [|.[] 

To illustrate the coincidence of transitions in bosonic and fermionic observables more 
precisely, we show in fig. |]the observables E± and (xx) an d also amp, all obtained in high 

1 The positions of the 1 st order transitions to the right and below the point T in fig. [| have been determined 
in a similar way. 
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Figure 3: Energy E± (eq. (p.lOj) ), (xx) and amp as functions of k at j3 = 0.70 for 
arriQ = 0.04. The run at k = 0.283 contains phase flips. 



statistics runs (1600-9600 trajectories), plotted as functions of k at (3 — 0.70 for am = 0.04. 
Again, the positions of the transitions agree within their accuracy on a finite lattice. At 
K = 0.283 we have observed 2-3 phase flips in E± and (xx), and the measured values of the 
observalbes are thus to some extent accidental. 

In order to demonstrate that amp has a discontinuity also in the chiral limit we have 
linearly extrapolated amp, obtained in long runs, to amo = very close below and above 
the phase transition, as shown in fig. f|. The value of amp at (3 = 0.70 on the 6 3 16 lattice 
drops from ~ 1 in the Nambu phase to nearly in the Higgs phase. 

The agreement of the positions of discontinuities of bosonic and fermionic observables 
suggests that the positions of the Higgs and chiral phase transitions on the ET line coincide. 
We conjecture that this is so, but we have to admit that the data leave open a logical 
possibility that the fermion observables do not jump exactly to zero as k increases across the 
Higgs phase transition. The genuine chiral phase transition might take place at some larger 
k, with small but finite values of fermionic observables in between, hidden by finite lattice 
size, finite amo, etc., effects. Therefore the - very crucial - conjecture of the coincidence of 
both transitions should be always checked in future studies with larger resources. 



4 NE line 

The NE line is apparently of second order, as we have found no indications of metastability 
in the range < (3 < 0.60. This was to be expected for small (3 where the chiral transition is 
essentially that of the NJL model J7| . As f3 increases the chiral transition gets substantially 
steeper, as demonstrated in fig. ^. 

At j3 = 0, we have checked that our data agree within good accuracy with high precision 
data obtained recently p(J for the NJL model. This is not completely trivial, as we use a 
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Figure 4: Fermion mass as a function of amo at (3 = 0.70 for k = 0.28 and n = 0.29 on 
the 6 3 16 lattice with a linear extrapolation to amo = 0. 
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Figure 5: Comparison of the dependence of the fermion mass arriF on k for various (3 > 0. 
k c is the extrapolated critical point at arriQ = in infinite volume (see table 111) . The data 



for (3 = (crosses) are taken from ref. [2C]. 
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Figure 6: The observable $7 1 defined in eq. ( PH ). The curve is a power law fit described 
in the text. 



form of the Lee-Shrock transformation, which is different from that applied in [B(J (no scalar 
field is used there, and the model after the transformation is not gauge invariant). 

Our aim was to find out how far the critical behavior in the vicinity of the NE line is 
analogous to that at (3 = 0, i.e. to the NJL model, when (3 is increased. That this must be 
so for some finite interval in (3 around (3 = follows from the strong coupling expansion [|7| . 
Our strategy is to compare our data at various (3 > with the data obtained in ref. 0] at 
(3 = using the same analytic Ansatze. 

The position of the phase transition has been estimated on the 6 4 lattice using the 
microcanonical fermionic average method f24 [, allowing a calculation at amo = 0, but without 
an extrapolation to the infinite volume. The method is described in the appendix. The 
observable Q, eq. (|2.11|) , obtained by this method at fixed k = 0.80, is shown in fig. ||. The 
positions of the pseudocritical points (3 pc on 6 4 lattice have been determined by means of a 
power law fit of the form VL~ l {(3) = A{(3 — (3 pc ) J to the significantly nonzero values of Q^ 1 . 
The obtained values of 7 are consistent with 1 within the systematic error ±0.2, estimated 
by fits to different sets of points. The resulting (3 pc are indicated in fig. [I] by crosses. 

The position k c {(3) of the chiral phase transition at am = on the infinite volume 
lattice has been determined by means of the hybrid Monte Carlo algorithm at finite arriQ 
and on various lattice sizes. The extrapolation to the infinite volume and m = has been 
performed by means of a modified gap equation at (3 = 0, 0.55 and at K = 0.60 (3 circles 
most to the left in fig. P as explained below. The results of the extrapolation are listed in 
table 0. 



The modified gap equation has been used in ref. for analytic description of the data 
for the chiral condensate in the NJL model on different lattices and for different amo. It is 
defined as follows: ^ 

®ti = vY. r2^2^ —T~ k ( 4 - 1) 



sin 



with 



J = <$w(k,/3)(xx) + c 2 am 



(4.2) 



This equation is a generalization of the usual gap equation by introducing free parameters 
w(k,/3), c\ and c 2 independent of am and lattice size. The c\ and c 2 are expected to be 
nearly independent of k for k ~ k c . In the original gap equation for the NJL model their 
values are w(k, 0) = 8G(k), C\ = r(«) and c 2 = 1 (for the definition of G and r see eq. ( p.6|) 
and ( ^.7| )). For k ~ k c the variation of r(/c) is small at /3 — 0. To use the modified equation 
also for /3 > we neglect this presumably weak dependence and keep C\ constant. Up to 
this simplification and some redefinitions of the parameters the equation Q4.1|) is the same as 
eq. (3.4) of ref. [p0|1 . We can thus compare our data analysis at (3 > with that performed 



in |0| at (3 = 0. 

We have used eq. ( |4.1| ) as an implicit equation for w. For each value of (3 we choose ci 
and c 2 , and then for each k ~ k c , amo and lattice size we determine w from the value of 
(xx)- The choice of c\ and c 2 , independent of k, amo and lattice size, is made so that w is at 
fixed (3 independent of am and lattice size as much as possible. We minimize the expression 

(wi - W(k,Ci,c 2 )) — (4.3) 
with respect to variations of c\ and c 2 . Here i«j and <7j are, respectively, values and errors 



of w obtained by means of (pTT| ) from value and its statistical error at different data 
points i. w is the average of Wj for the same k, c\ and c 2 . 

The possibility of choosing c\ and c 2 in such a way that w depends only on k means that 
the modified gap equation is applicable at this (3. The results for w(k) are shown in fig. [7]a,b. 
If the modified gap equation works, it can be used for an extrapolation to chiral limit and 
infinite volume, and the critical value w c of w, below which (xx) vanishes (see fig. 0), is 
obtained. The critical point K c (/3) at the given (3 is then determined from w(k c ) = w c . 

As seen from fig. [7|a and 03, the modified gap equation describes in the above sense 
the data very well at (3 — and (3 = 0.55 for various amo and different lattice sizes. The 
agreement at /3 — 0.55 (fig. H) is nearly as good as at (3 = 0, cf. ref. |p0[] . The constants c\ 
and c 2 have nearly the same values for these two (3 values. Therefore this method can be 
used also for fixed k and varying (3 ~ (3 C . We have done this at k = 0.60. The lowest curve 
in fig. | describes the extrapolation to infinite volume and amo — 0. In this way the position 
n c (f3) of the chiral transition has been obtained for 3 points on the NE line. 

In this way one can obtain also pseudocritical points on lattices of finite size, again by 
determining values of w for which (xx) vanishes. The corresponding (3 pc on the 6 4 lattice 
are consistent with those obtained by the microcanonical fermion average method. 

Another way to investigate the scaling properties of the NJL model in PUI was to plot 
(xx) 2 as a function of am /(xx) (Fisher plot). This is motivated by the general scaling 
relation 

am = 6i(« - « C )(XX>^"^ + Hxx) 6 , (4-4) 

where the mean field values of the exponents, 6 = 3 and (3 = ^, are assumed. Provided these 
values are correct, the data for a fixed k on a given lattice should form straight lines. As 
seen in fig. pla, for f3 = 0.55 the mean field exponents describe the data quite well. We again 



find a large similarity between our data at f3 = 0.55 and the data at (3 = in ref. p0| . This 
includes also the observed finite size effect. 
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Figure 7: The parameter w in the modified gap equation as a function of k for (a) = 
(NJL model), (b) (3 = 0.55 and (c) = 0.64. For = and = 0.55 w is independent 
of lattice sizes and bare mass and the modified gap equation describes the data well. For 
= 0.64, which is close to the point E, this is no longer so. The curves in (a) and (b) are 
quadratic function fits to the u>-values. The data in figure (a) are taken from ref. 
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Figure 8: Comparison of the modified gap equation and the data at j3 = 0.55. The 
parameter w depends on k as described by the line in fig. ^b, whereas the parameters c\ 
and C2 are fixed at values c\ = 0.74 and C2 = 1.1. The extrapolation to infinite volume and 
arriQ = is represented by the dashed line. 



5 Vicinity of the point E 

It is notoriously difficult to localize a tricitical point. Our best current estimate is /3g = 
0.64+q;q4 and k e = 0.305^q;^ 5 (point E in fig. [1]). It is based on the extrapolation in (3 of 
latent heat to zero and on the localization of minimum of scalar mass arris as described 
below. Finite size effects are only roughly taken into account. 

We have tried to apply at (3 = 0.64 the same method of analysis to the scaling behavior 
as at (3 = and (3 = 0.55. It turns out that both above methods do not work anymore. 
The modified gap equation does not describe the data if w should depend only on k and not 
on arrio and lattice size. This is seen in figure |7|c. Here we show the values of w obtained 
for such a set of parameters c\ and C2, for which the dependence of w on arriQ and lattice 
size is minimal. Also the test for mean field values of exponents j3 and 5 by means of the 
Fisher plot fails. In fig. |b, the lines connecting data at the same k and lattice size differ 
significantly from straight lines. 

We have found another difference between the behavior of the model close to the chiral 
phase transition around the point E and at smaller (3. The scalar boson mass arris gets 
small in the vicinity of E, whereas it remains of the order of cutoff or larger at the chiral 



transition for smaller (3. In fig. [10], both boson masses arris and amy are shown at (3 = 0.55 
and (3 = 0.64 for arriQ = 0.04. A dip in arris at /3 — 0.64 around k = 0.315 is clearly 
seen. It gets deeper when the lattice size is increased. Similar dips have also been found for 
other values of am and lattice sizes for various (3 ~ /3e, whereas at (3 = 0.55 or smaller (3 
the mass arris does not decrease below 1 neither on small lattices nor on large ones. This 
difference, obvious in the quenched case, is thus present also with dynamical fermions. If a 
new correlation length diverges at some critical point and not at other ones, the most natural 



Fisher plot 8 0.55 



20 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 V |- I 1 1 1 




18 


- 






16 






6 3 16: 


14 






X:/c = 0.37 
o:/c = 0.38 


12 






□ :/c=0.39 


10 


® 




*:/c=0.4 
n:/c = 0.41 


08 






8 3 24: 


06 






O:/c = 0.37 


04 






A;/c = 0.38 
*:/c = 0.39 


02 


i i i A i i i i i i i i 


ptttI ttttIttti 














.05 .10 .15 .20 .25 



am /<xx> 



Fisher plot ,3 = 0.64 



(b) 



A 
X 
iX 
V 



.20 
.18 
.16 
.14 
.12 
.10 
.08 
.06 
.04 
.02 





.05 



.10 



.15 



.20 



.25 



am 



0/ 



6 3 16: 
X:/c = 0.31 
o : /c = 0.313 
□ :/c = 0.32 
C:/c = 0.33 

8 3 24: 
n:/c = 0.30 
O:/c = 0.31 
a : /c = 0.313 
*:/c = 0.315 
^:/c = 0.32 



Figure 9: (xx) 2 as a function of amo/(xx) f° r fixed coupling parameters and lattice size 
(Fisher plots) at (a) (3 = 0.55 and (b) (3 = 0.64 ~ (3%. The straight lines in (a) indicate 
according to eq. ( ^4.4|) a mean field behavior. 
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Figure 10: Vector and scalar boson masses, amy and arris, as functions of k at (a) (3 = 0.55 
and (b) [3 = 0.64 ~ /3e- To estimate the finite size effects two different lattice sizes (6 3 16 
and 8 3 24) are shown. 
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Figure 11: Meson states observable in the vicinity of the chiral phase transition. The 



states are denoted according to ref. [29|. 



explanation is that the scaling behavior is different. 

These observations indicate that the scaling behavior in the vicinity of the point E is 
different from that at smaller f3. This tentative conclusion is supported by the experience 
with tricritical points in various statistical mechanics models f31fl . The tricritical exponents 
have usually quite different values from those describing the universality class of the adjacent 
critical manifolds. We also note that the chiral and magnetic phase transitions, occurring in 
the xU<fii model at small and large (3, respectively, meet at the point E. Thus phenomena 
quite different from those known in both limit cases might occur here. 



6 Summary and outlook 

From the study of the behavior of several observables in the vicinity of the point E, as 
well as from some more general considerations about tricritical points, we conclude that the 
tricritical point E in the phase diagramm of the xU<p4 model is an interesting candidate for 
an investigation of the continuum limit of this model. The significant differences from the 
NJL-type model suggest that a continuum theory with massive fermions might arise which 



does not belong to the universality class of the Yukawa models in the sense of ref. |2T . 

To test this conjecture a laborious investigation of the scaling behavior of the model 
around the tricritical point E will be necessary, taking more observables into account. A 
rich spectrum is expected in the Nambu phase of the model, qualitatively similar to that one 
would guess in the case of QCD with one fermionic and one scalar quark (this analogy has 
been discussed in ref. Some glimpse at the states not too difficult to find is provided 
by figs. [TO] and 11]. In fig. |TTJ we show some meson masses. It has been checked in []lj| that 



for the Goldstone boson mass am^ the relationship am n oc A /amo, expected in the Nambu 
phase, holds. In this figure one can identify also states which are analogies of the p, a and 



A\ mesons in QCD. Such states might accompany the massive fermion F and the scalar 
S also in the continuum limit of the xU<fii model constructed at the tricritical point E. If 
so, the spectrum would be rather different from that of the Yukawa models, supporting the 
conjecture of a new universality class. 
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A Microcanonical fermionic average method extended 
to field theories with scalars 

To perform the calculation of fi _1 , eq. ( |2.11| ), in the chiral limit the microcanonical fermionic 
average method developed for the gauge- fermion systems |22],^3|] has been extended to the 
fermion-gauge-scalar models Here we just outline the strategy. 

The essential idea of the algorithm is the computation of the full effective action 

S eS {E,m ,N f ,p,K) = -IuM{E,k) -6V(3E + S^{E,m ,N f ,K) (A.l) 

as a function of the pure gauge energy E and other bare parameters like m , Nf, [3 and k, 
where 

M(E, /t) = J [dU}[d<f)}[dtf}5[Ep{U) - E] exp[8VKE L (U, 0)] (A.2) 
is the density of states and 

/ N f 

S^(E,m ,N f ,K)\ = -ln/[detA(f/)]- 



N r 



^ J[dU][d<j)}[d^]6[E P (U)-E} [detA(U)]^exp[8VKE L (U)} 

is the effective fermionic action. M(E, k) in (|A.2|) can be easily computed by quenched 
simulation. The main effort to be paid is to calculate S[ s (E,m , Nf, n) through ( |A.3| ) by 
means of microcanonical and Monte Carlo simulations For Nf = 4, (det A(U))e is the 
fermionic determinant averaged over the configurations with the probability distribution 



5[E P (U)-E}ex P [8VnE L (UA)] f , 

M{E~Kj • (A - 4) 

From (|A.3|) , one sees that the effective fermionic action S^ s does not depend on (3 and its 
dependence on m and Nf can be easily obtained, once the eigenvalues of the massless 
fermionic matrix for the decorrelated configurations at fixed E have been calculated by the 
Lanczos algorithm. Once S e ff(E,mo, Nf, (3, k) is known, the thermodynamical quantities 
such as (E P ), (E L ), {tptp) and can be obtained by deriving the partition function 



Z= I dEexp[-S eS (E,m ,Nf,P,K)} 



(A.5) 



or by saddle point analysis . 



In [pi], on 6 4 lattice and at finite arriQ, comparisons with the results from the Hybrid 
Monte Carlo (HMC) method for a fermion-gauge-scalar model were made, and the HMC data 
were well reproduced. The advantage of the algorithm, in addition to its lower computational 
costs in searching the parameter space (am Q ,/3) compared to the conventional algorithms, 
is the accessibility of the chiral limit. Therefore the algorithm is very suitable for the phase 
structure analysis. 

While having demonstrated the advantages of the algorithm on the 6 4 lattice, we would 
like to mention some open questions. For large lattice volumes, it is time and memory 
consuming to diagonalize the fermionic matrix. For example, on the 8 4 lattice, the amount 
of data required grows so fast that the algorithm has not been well tested. Although the 
Lanczos algorithm is parallelizable ||23|| , it is not easily vectorizable. It is also not easy to 
evaluate observables other than the thermodynamical quantities. These problems are under 
further investigation. 
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